Draft version November 1,2011 

Preprint typeset using LATgX style emulateapj v. 26/01/00 



O 

(N 

o 
O 

cn 
< 

Oh- 
I 

o 

(N 
> 

(N 



SUBSTRUCTURE IN THE STELLAR HALOS OF THE AQUARIUS SIMULATIONS 

Amina Helmi\ a. p. Cooper^'^, S.D.M. White^, S. Cole^, C.S. FRENK^ J.F. Navarro'* 

Draft version November 1, 2011 

ABSTRACT 

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by 
the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius 
Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams 
akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a 
characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS 
of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated 
for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely 
to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ^ 10% 
of the Milky Way halo stars have formed in-situ. 

Subject headings: galaxies: halos, structure, formation, - Galaxy: halo, kinematics and dynamics, structure 



1. INTRODUCTION 

Stellar halos are repositories of merger debris and hence, 
despite their low luminosity, are central to efforts to unravel 
the accretion history of galaxies jSearle & Zinn|[T978l) . Wide- 
field photometric surveys have discove red a plethora of sub- 
structures in the halo of our Galaxy (e.g . iMaiewski et alj|2003t 
iBelokurov et al."2006'; 'Juric et al.'^OOS*). Some examples are 
the Orphan Stream (Belokurov et al. 2007a) and the Sagittarius 
dwarf (Ibata et al. 1994) and its extended tidal tails dlvezic et al.l 
12000' .Yannv et al., 2000) . Other conspicuous stru ctures, like 
the Hercules- Aquila cloud (IBelokurov et alj|2b07bh . the Virgo 
overdensity and the Virgo Stellar Structure, do not have the 
elongated appearance typical of tidal features and hence their 
nature is less obvious. Moreover, some of these structures 
strongly overlap on t he sky but their relationship is not al- 
ways clear (Mai-tmez- Delgadoetal.l 120071: iDuffau et al] 120061: 
iNewberg et al. 2007). 

The amount of spati al substructure pr esent in the outer halo 
has been quantified bv lBell et al.l ( l2008l see also de Jong et al. 
2010). These authors found RMS fluctuations of order 30^0% 
with respect to a smooth halo model. Taken at face value, this 
implies that at le ast this fr action of the stellar halo has been ac- 
creted. Starkenburg et al.l (|2009 ) determined that at least 10% 
of the outer halo must have been accreted, based on the amount 
of clustering present in the Spaghetti Survey (which also in- 
cluded radial velocity information). Roughly half of this sur- 
vey's pencil beams pointed to portions of the sky unexplored 
by SDSS; intriguingly, these authors found no evidence of sub- 
structure in those regions, suggesting that the distribution of 
outer halo substructure may be very anisotropic. 

Several studies have examined the properties of stellar ha- 
los in ACDM simulations of s tructure formation. For exam- 
ple, iBuUock & JohnstonI (|2005|) followed the accretion of satel- 
lites onto a smooth, time-dependent gravitational potential re- 
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sembling a Milky Way galaxy. Their simulations produce 
very lumpy stellar halos, with extended tidal stream s covering 
large portions of the sky. iDe Lucia & Heknil (l2008h used high- 
resolution cosmological N-body simulations of the formation of 
a Milky-Way size dark halo in combination with a semi-analytic 
model of galaxy formation. Although their model reproduces 
the global properties of the Galactic stellar halo, the resolution 
of the simulations was too low to show much substructure in 
the form of tidal tails. 

Recentlv. lCooper et al.l (1201 Ol hereafter CIO) have combined 
full cosmological dark matter simulations from the Aquarius 
project with the Durham semi-analytic galaxy formation model, 
GALFORM, in order to follow the assembly of 'accreted' stellar 
halos. The extremely high resolution of the Aquarius halos per- 
mits a wealth of substructure to be studied in the Cooper et al.l 
models. This Letter will focus on the properties of this sub- 
structure as may be revealed by ongoing (e.g. SDSS) and fu- 
ture wide-field photometric surveys, such as Skymapper, Pan- 
STARRS or LSST 

2. BRIEF DESCRIPTION OF THE MODEL 

The Aquarius project consists of simulations of the formation 
of several dark matter halos with masses comparable to that of 
the Milky Way in a ACDM univers^. The Aquarius halos (la- 
belled Aq-A to Aq-E) are resolved wit h more than 10^ parti - 
cles, and their properties are described in ISpringel et ai](l2008l) . 
To follow the evolution of the baryonic component, the dynam- 
ical information obtained from the simulations is c oupled to the 
Durham semi-analytic galaxy formation model dBower et al.l 
120061) . In essence, CIO identify at each output the halos that 
host galaxies according to the semi-analytic model, and tag the 
1 % most bound particles at that time with "stellar properties" 
such as masses, luminosities, ages and metallicities. This al- 
lows us to follow the assembly and dynamics of the accreted 
component of galaxies. 

Note, however, that we do not follow stellar populations 
formed in-situ (for example, a disc component) because our 
dynamical simulations are purely collisionless (hydrodynam- 
ics and star formation are treated semi-analytically). A massive 
disc may influence satellite orbits and their debris, and also pro- 

^^Here we focus on the analysis of five of the six Aquarius halos, because 
the remaining object's entire stellar halo was built in a merger at 2 ~ 0.3 and 
is thus unlikely to be representative of the Milky Way. 
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vide a reservoir from which stars can be scattered into the halo, 
but these effects are absent in our models. 

The tagging scheme leads to the identification of more than 
4 X 10^ tracer dark matter particles for Aq-A and up to a 
maximum of 6 x 10^ for Aq-B and Aq-D. Since the stellar 
mass assigned to each of these dark matter particles depends 
on the mass of stars formed in the parent galaxy with time, 
they do not all carry the same weight (for more details, see 
CIO). Therefore, and for ease of comparison to observations, 
we have resampled the tagged particles as follows. For each 
tagged dark matter particle we find the 32 neighbors ranked as 
nearest simultaneously in position and in velocity at the present 
time amongst the class of tagged particles that share the same 
progenitor We then measure the principal axes of the spatial 
distribution of these neighboring particles. This is used to de- 
fine the characteristics of a multivariate Gaussian from which 
we draw positions for re-sampled stars. We then down-sample 
this new set of "stellar particles" to represent main sequence 
turn-o ff stars (assuming 5 Mq per MSTO, as used in lBell et"aLl 
|2008|) or red giant branch stars (we assume 1 R GB per 8 MSTO 
stars, as in the models of Mari go et al.ll2"008h . This results in 
~ 5.5 X 10'' MSTO stars between 1 kpc and 50 kpc from the 
halo center for Aq-A, and up to 1.83 x 10^ MSTO stars for 
Aq-D (which is the brightest of our stellar halos). 

3. RESULTS 

In this Letter we focus on the distribution of halo stars on 
the sky for direct comparison to photometric surveys. We will 
show that panoramic views resembling the "Field of Streams" 
are common in our simulations, and we will provide a new 
theoretical perspective on the nature of the overdensities dis- 
covered by SDSS. At the end, we quantify the amount of sub- 
structu re present in our stellar halos and compare it to the re- 
sults of lBell et alj (1200 8^. In future papers we will focus on the 
characterization of substructure (width, age, relation to progen- 
itor system, etc.), and on the amount of kinematic substructure 
present in the Solar neighborhood (especially for comparison 
with RAVE, and Gaia data in the future; Gomez et al. in prep.). 

3.1. Distribution on the sky 

Figure [T] shows the distribution in the sky of the RGB stars 
present in the stellar halo of Aq-A. The different panels in 
Fig. [U correspond to stars located at increasing distance from 
the Sun, which is assumed to be at (-8, 0, 0) kpc with respect 
to the halo center (for this strongly prolate halo, the major axis 
is very nearly aligned with the z-axis of the simulation refer- 
ence frame). Particles that have the same color represent stars 
formed in the same parent galaxy. 

This figure shows that the distribution of stars in the accreted 
halo is very smooth in the inner 10 kpc. Substructure becomes 
apparent at ~ 20 kpc and dominates the halo beyond ^30 kpc. 
The substructures are often diffuse, particularly at small dis- 
tances from the galactic center. This is because their progenitor 
satellites are relatively massive. For example, the most promi- 
nent streams in Fig. [T] are those in magenta (visible at all dis- 
tances), green (dominant in the very center), blue-gray (which 
we describe below as a Sgr analogue) and light green (promi- 
nent beyond 30 kpc). These conti'ibute 1.3 x 10^, 1.4 x 10^, 
2.2 X 10'' and 4.5 x 10'' M© respectively, i.e. all together they 
add up to 85% of the stellar mass in the halo of Aq-A. 

The distribution of substructure on the sky is anisotropic, and 
appears to be preferentially found along a "ring". In this re- 
gion, a system of streams very similar to those of the Sagittarius 



dwarf galaxy is apparent (shown in blue-gray). Even the central 
regions of the stellar halo of Aq-A are not isotropic on the sky 
- a bar-like feature is visible towards the galactic center The 
central region of this stellar halo i s triaxial, which may also be 
true in the case of the Milky Way dNewberg et al.ll200'7l) . How- 
ever, it is likely that the degree of triaxiality would decrease if 
a massive stellar disk were included in the simulation (see e.g . 
'Dubinski 1994; Debattis ta et"ani2008l; iKazantzidis et alj|2004 
,2010; Tissera et al. 2009) . 

Figure 12] shows the sky distribution of stars located between 
10 and 30 kpc (top) and 30 to 50 kpc (bottom) for the remaining 
halos. Substructure is apparent in all cases, but there are large 
variations in its coherence, surface brightness and distribution 
across the sky. As in the case of Aq-A , the most significant 
overdensities are found at 10-30 kpc, although abundant sub- 
structure is also present at larger distances. 

The object found in Aq-A that resembles Sagittarius and its 
streams became a satellite at redshift 1.74 (9.7 Gyr ago), and 
had a total mass of 2.9 x 10^° and a stellar mass of 2.7 x 
10*^ M0 at the time of infall. A bound core survived until the 
present day with approximately 25% of the initial stellar mass. 
Fig.[3]shows that its streams cover a similar area of the sky and 
have distances comparable to those of the Sagittarius streams in 
the Milky Way halo. 

We also find structures that resemble the Orphan stream. Al- 
though they are not common at small radii, very low surface 
brightness thin streams are present, and can be found as close 
as 10 kpc from the halo center. Typically they have elongated 
orbits, and their characteristic apocentric distances range from 
30 kpc up to ^ 130 kpc. This explains why they have remained 
coherent despite the triaxial shape of the halo (and the grain- 
iness of the potential). Such streams originate in low mass 
galaxies. For example, the progenitor of the thin stream shown 
in Fig. [3]had a total mass of 1.2 x 10* Mq, a stellar mass of 
10^ Mq and was accreted into Aq-A at z = 5, i.e. ~ 12.5 Gyr 
ago. 

3.2. The "Field of Streams" 

The previous figures show that in our simulations, streams 
from different progenitors frequently overlap on the sky. This 
is due mostly to the correlated infall pattern of the progeni- 
tor satellites from which these streams originate. For exam- 
ple, Aq-A remains embedded in the same large-scale coherent 
filament for ^ 10 Gyr before the present day, and many of its 




Fig. 3. — Sky distribution of RGB stars in Aq-A from tlie "Sagittarius" 
(blue) and "Orplian" (orange) streams analogues. Only those stars located at 
distances smaller than 80 kpc from the "Sun" are shown. 
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Fig. 1 . — Distribution of "field" RGB stars on tlie sky at various distances from tlie "Sun" for the stellar halo of Aq-A. The different colors correspond to stars 
originating in different progenitors. The total number of progenitors is 163, out of which 53 contribute more than 10* Mq in stars. There are 35 different progenitors 
contributing at least 50 RGB stars in the innermost bin, and 48 at distances between 10 and 30 kpc. This number drops down to 10 in the most distant bin considered 
here. Thus the Aq-A stellar halo is most diverse between 10 and 30 kpc. 



satellites have formed in this structure fcf lLibeskind et al.l2005L 
also Lovell et al. 2010, Vera-Ciro et al. in prep.). 

This is exemplified in Fig. H] which plots the distribution of 
MSTO stars (generated following our simple prescription) for 
Aq-A in a thin distance slice through a region similar to the 
SPSS footprint known as the "Field of Streams" (iBelokurov et al.l 
l2006h . Here different colors indicate different progenitor sys- 
tems. The characteristics of this Aquarius "Field of Streams" 
are similar to those observed in SDSS. In particular, we see 
streams of stars that follow si milar paths on the sky, resembling 
the bifurcation discovered bv lFellhauer et al.l (l2006l) . as well as 
various broad overdensities. Note that some of the bifurcations 
that are apparent do not necessarily arise from the overlap of 
streams on the same orbit with different orbital phase, but in- 
stead correspond to the overlap of streams of different origin. 
This implies that measurements of position and distance alone 
may be insufficient to associate overdensities in nearby regions 




Fig. 4. — Distribution of MSTO stars in the Aquarius "Field of Streams" 
colored according to provenance. The region of the sky shown here has a sim- 
ilar extent to the SDSS footprint, but the model "stars" are located in a thin 
slice of 1 kpc width at a distance of ~ 35 kpc. Streams from different pro- 
genitors overlap on the sky because of the correlated infall directions and also 
due to group infall (both characteristic of ACDM). This constitutes a plausible 
(alternative) explanation of the observed bifurcation of the Sgr stream. 



of the sky with a single parent object; therefore some caution is 
required when such associations are used to constrain the shape 
of the underlying gravitational potential. 

3.3. Quantifying the amount of substructure 

Following iBell et al.l (l2008h . we have selected stars from 
SDSS-DR7 with 140° < a < 220° and 0° < ^ < 60°, 

0. 2 < [g-r) < 0.4 (characteristic of the halo MSTO), and with 
apparent magnitude 18.5 <r< 22.5. If an absolute magnitude 
Mr ^ 4.5 is assumed, these stars would probe a distance range 
of 7 up to 35 kpc. In practice, this simple 'tomography' of the 
halo is subject to various uncertainties. 

For example, using data from SDSS Stripe 82, lJuric et all 
(2008) estimated that at the bright end, 5% of the point- like 
sources with 0.2 < {g-r) < 0.3 are QSOs. At r 22, this 
fraction increases to 34%, which extrapolated to r ^ 22.5, im- 
plies a non-negligible contamination of 50%. Another com- 
plication arises from photometric errors. Although these are 
relatively small, the color error at the faint end is ag^r ~ 0.2, 

1. e. co mparable to the width of the MSTO color selection box 
dlvezic et al..,2008i) . Therefore, a large number of lower main 
sequence stars are scattered into the MSTO region, leading to 
an incorrect assignment of absolute magnitudes. We consider 
both of these effects in our simulations. 

We measure the RMS in the sky density of sources identi- 
fied as halo MSTO stars in the SDSS as a function of apparent 
magnitude, in cells of 2x2 deg^. To compare with the simula- 
tions, we derive the apparent magnitude of each model MSTO 
stellar particle as m,. = 5 logd — 5 + Mr + AMr- Here d is 
the distance to the observer, and Mr is a random variable fol- 
lowing a Gaussian distribution with (Mr) = 4.5 and a = 0.9 
(to mimic the spread in the absolute magnitude of the MSTO). 
The term AMr accounts for the uncertainty due the color error 
(Jg-r, and is computed as AMr — <i^Ir/d{g — r) x Ug-r x u, 
where dMr/d{g — r) is the slope of the main seq uence track 
(derived from the photometric parallax relations by 'Juric et all 
I2OO81) . and u is a normally distributed random variable with 
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Fig. 2. — Sky distribution of RGB stars located at "lieliocentric" distances [10, 30] kpc (top) and [30,50] kpc (bottom) for the set of Aquarius stellar halos. 
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Fig. 5. — RMS to mean density of MSTO "stars" inaregion of the sky sim- 
ilar in location and extent to the SDSS footprint as a function of apparent mag- 
nitude for our accreted Aq-halos. Top panel: The dashed curves indicate the 
raw measurements, while the solid curves take into account photometric errors 
and QSO contamination. Bottom panel: A smooth component with 10% of the 
total stellar halo mass has been added to each Aq-halo. The measurements for 
eight observers located along the "solar" circle are given by the dotted curves, 
while their median value is shown in solid. In both panels, the black asterisks 
correspond to the values estimated from SDSS data for the Milky Way. 



zero mean and unity dispersion. 

In our simulations, we place eight observers along a circle of 
8 kpc radius from the halo center, in a plane perpendicular to 
the minor axis of the dark matter halo (an approximation to the 
Solar circle). We select those MSTO "stars" located in the same 
region of the sky and with the same apparent magnitude range 
as the SDSS sample, and we measure the (Poisson-corrected) 
RMS in the projected surface density of these stars for each of 
our halos. Finally, we add the expected contamination by QSOs 
as described above. In this way, we obtain the fractional RMS 
which may be compared to the measurement from SDSS. 

In the top panel of Fig. |5] the solid curves denote the me- 
dian (of each of the 8 independent observers) for each halo. 
The dashed curves correspond to the raw RMS obtained with- 
out taking into account the effect of photometric errors and 
QSO contamination. This figure shows that our halos tend to 
have systematically larger fractional RMS values than found in 
the SDSS data. The bottom panel of Fig. |5] shows the RMS 
statistic when we include the contribution of a smooth compo- 
nent. We plot both the median RMS for each halo (solid) as 
well as that measured by each of the 8 observers (dotted). Here 
we have assumed 10% of the total stellar halo mass to be dis- 
tributed following a Hernquist profile with scale radius of 1.25 
kpc, although our results are not strongly dependent on the ex- 
act value of this parameter (note that at r '--^ 19, its contribution 
is ~ 30%). This component would represent stars formed in the 
Galaxy itself, which, by definition, are absent in the CIO mod- 
els of accreted halos. The main effect is to lower the contrast 
between cells with and without overdensities. The comparison 
to the SDSS results improves noticeably, indicating that there 
is indeed room for an in-situ compon ent to be present in the 
stellar halo of the IVIilky Way (see also lZolotov et al.ll2009l) . 

4. SUMMARY 

We have discussed the substructure present in the stellar ha- 
los of Cooper et al. (2010) which were formed by the disrup- 
tion of accreted satellites within the IVIilky Way mass ACDM 
halos of the Aquarius project. Diffuse features such as the 
Pisces and Virgo overdensities and the Hercules-Aquila cloud 
are relatively common, as are narrow, faint "orphan" streams. 
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One of our halos even contains a structure similar to the Sagit- 
tarius Stream. We have found that substructures are not dis- 
tributed isotropically on the sky. This is likely to be related to 
the correlated infall of satellites accreted along large-scale fil- 
aments, a characteristic feature of the cold dark matter m odel 
(iLibeskind et al.ll2005l: lUfe Helmill2008l: iLux et al.ll20Tol) . As 
a consequence, chance aUgnments between streams from dif- 
ferent progenitors may occur (as in our "Sagittarius" example), 
which may be easily misinterpreted as streams from the same 
object with different orbital phase. We find that the fraction of 
substructure in our models is somewhat higher than is estimated 
for the Milky Way. The addition of a smooth spheroidal compo- 
nent that contains 10% of the total stellar halo mass is required 
to match the simulations the SDSS data. This might suggest 
that a fraction of the halo stars formed in-situ. However, the 
discrepancy could also, however, be the result of limitations in 
our modeling techniques or uncertainties in the interpretation 
of the SDSS data. 
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